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^vj ■ The addition of non-renormalizable terms involving the Higgs fields to the MSSM ameliorates the little 

hierarchy problem of the MSSM. For neutralino dark matter, new regions for which the relic abundance 
of the LSP is consistent with WMAP (as the bulk region and the stop coannihilation region) are now 
permitted. In this framework, we analyze in detail first the direct dark matter detection prospects 
in a XENON-like experiment; then we study the capability of detecting gamma-rays produced in the 
Q , i annihilation of neutralino LSPs in the FERMI satellite. 
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D . 1 Introduction 



The smallness of the quartic Higgs coupling in the framework of the minimal supersymmetric standard 

model (MSSM) poses a problem. The tree level bound on the Higgs mass is violated, and large enough 
J> ' loop corrections to satisfy the lower bound on the Higgs mass suggest that the stop sector has rather 

^^O peculiar features: at least one of the stop mass eigenstates should be rather heavy and/or left-right-stop 

mixing should be substantial^. 
f»o ■ The situation is different if the quartic Higgs couplings are affected by New Physics. If the New 

. * | Physics appears at an energy scale that is somewhat higher than the electroweak breaking scale, then its 

effects can be parametrized by non-renormalizable (NR) terms. The leading NR terms that modify the 
^^ , quartic couplings are^: 

Wbmssm = J^( H uH d ) 2 + jjZ(H u H d ) 2 , (1) 

where Z = 9 2 m susy is a SUSY-breaking spurion. The first term in equation [T] is supersymmetric, while 
the second breaks supersymmetry. In the scalar potential, the following quartic terms are generated: 

2e 1 H u H d {HlH u + H d H d ) + e 2 (H u H d ) 2 , (2) 

where 

61 s -3JT' £2 " —■ (3) 

Let us note that, in some regions of the parameter space, these operators Beyond the MSSM (BMSSM) 
may destabilize the scalar potential - . 

One of the attractive features of the MSSM is the fact that the lightest SUSY particle (LSP), usually 
the lightest neutralino, is a natural candidate for being the dark matter (DM) particle. The effects of the 
NR operators are potentially important in the determination of regions fulfilling the WMAP DM relic 
density constraint . On the one hand, these operators give rise to a new interaction Lagrangian which 
contributes to the vertex of two higgsinos and one or two Higgs bosons. This effect is relevant when the 
neutralino LSP has a significant component of higgsinos. On the other hand, the uplift of the lightest 
Higgs mass could reopen regions giving rise to a relic density in agreement wit the WMAP results, but 
which were ruled out by the boun ds ove r the Higgs mass. The effect of these operators on the relic density 
was studied in detail in references I I I In addition, the con sequences for baryogenesis and electric dipole 
moments in the BMSSM has been studied in references'^^. We shall evaluate the detection perspectives 
for two different detection modes, namely direct detection in a XENON-like experiment and gamma- ray 
detection from DM annihilations in the galactic center for the FERMI mission. 



2 The model 

The BMSSM framework, if relevant to the little hierarchy problem that arises from the lower bound on 
the Higgs mass, assumes a New Physics scale at a few TeV. Since the new degrees of freedom at this 
scale arc not specified, the effect of the new threshold on the running of parameters from a much higher 
scale cannot be rigorously taken into account. It therefore only makes sense to study the BMSSM effects 
in a framework specified at low energy. Within this framework, we calculate the DM relic density, and 
the direct and indirect detectionprospects in the presence of the new ti couplings. We used a modified 
version of the code micrOMEGAs *■, where we implemented the BMSSM couplings, in order to calculate 
the relic density as well as the cross-sections and decay channels relevant for DM detection. The leading 
e^-induced corrections to the spectrum, were implemented using the code SuSpectM^. 

The first scenario considered contains correlated stop and slepton masses, as the mSUGRA frame- 
work. In this case, the ncutralino LSP is an almost pure bino-like state; the 'bulk region' is highly 
constrained due to the experimental lower bound on the Higgs mass. In general, in order to fulfill such a 
constraint either heavy or mixed stops are required. Additionally to the ordinary mSUGRA parameters 
we have two extra BMSSM parameters: ei and e 2 . Let us emphasize again that one should not think 
about this scenario as coming from an extended mSUGRA model, since the effects of the BMSSM physics 
at the few TeV scale on the running cannot be taken into account. The upper panels of figure [1] show the 
regions in the [too, TO1/2] plane in which the WMAP constraint is fulfilled (red lines). We fixed Aq = 
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Figure 1: Regions in the [mO, mi/2] plane (upper panels) or [Mi, p] plane (lower panels) in which the WMAP constraint 
is fulfilled (red lines). The blue region is excluded by the null searches for charginos at LEP. The orange areas in the upper 
(lower) panels are excluded because the stau (stop) is the LSP. The dash-dotted black curves in the up panel are contour 
lines for the Higgs mass with values in GeV as indicated. Above the dotted purple line of the left of this line, the electroweak 

vacuum is metastable. 



GeV, fi > 0, tan/3 = 3, e 2 — and ei = 0, —0.1. The blue region is excluded by the null searches 
for charginos at LEP. The orange area is excluded because the stau is the LSP. The dash-dotted black 
curves are contour lines for the Higgs mass. Concerning the upper panels, a generic point in the former 
parameter space usually gives rise to a too large relic density, in conflict with the WMAP measurements. 
However, for moderate mi / 2 and low too values, there is a region where the LSP is almost degenerate in 
mass with the lightest stau, enhancing the co-annihilation cross-section X? — f. Another region giving 
rise to relic density in agreement with the WMAP measurements appears for TOi/ 2 ~ 120 GeV (top left 
panel). This is the '/i-pole' and the 'Z-pole' region in which to/j ~ mz ~ 2to x o, and the s-channel Higgs 
and Z boson exchange is nearly resonant, allowing the neutralinos to annihilate efficiently. Let us note 



that this region is already excluded by LEP measurements. For negative enough ei values (upper right 
panel), the uplift of the Higgs mass generates a splitting among the '/i-pole' and the 'Z-pole' regions, 
with the former now evading LEP constraints. 

The second model correspond to a pure low energy scenario giving rise to light unmixed stops. In 
addition to the BMSSM e, parameters, we consider the following set of parameters: 

Mi, \i, tan/3, X t , mjj, m Q , mj, m A , (4) 

where m i is a common mass for the sleptons, the first and second generation squarks, and bn. We further 
use Mi = | tan 2 6\y M^. To demonstrate our main points, we fix the values of all but two parameters as 
follows: ei = or -0.1, e 2 = 0, tan^ = 3, X t = 0, m v = 210 GeV, m Q = 400 GeV, mj = m A = 500 
GeV. This scenario gives rise to relatively light stops, namely m^ < 150 GeV and 370 GeV < m^ < 400 
GeV. We scan over the remaining two parameters, Mi and fi. 

The lower panels of figure [1] show the regions in the [Mi , /j] plane in which the WMAP constraint 
is fulfilled (red lines). The first region is the Z- and /i-poles in which the LSP is rather light, and the 
s-channel Z or h exchange is nearly resonant, allowing the neutralinos to annihilate efficiently. There is 
also a 'mixed region' in which the LSP is a higgsino-bino mixture, Mi ~ /i, which enhances (but not too 
much) its annihilation cross-sections into final states containing gauge and/or Higgs bosons. Finally the 
'stop co-annihilation' region, in which the LSP is almost degenerate in mass with the lightest stop. Such 
a scenario leads to an enhanced annihilation of sparticles since the Xi — *i co-annihilation cross-section 
is much larger than that of the LSP. 

Let us first consider the case where ei = e 2 = (left lower panel). In any case, for the Higgs mass 
values obtained here, m^ ~ 85 GeV, this region is already excluded by the negative searches for chargino 
pairs at LEP. The main difference when considering the ei = —0.1 case (lower right panel) comes from 
the important enhancement of the Higgs mass due to the presence of the BMSSM operators. In this 
case it is possible to disentangle the Z and the h peaks, since the Higgs-related peak moves to higher 
Mi values, due to the increase of the Higgs mass: mh = 122 GeV. Let us note that the latter peak is no 
longer excluded by chargino searches. 

3 Direct detection of dark matter 

DM direct detection experiments measure the number TV of elastic collisions between WIMPs and target 
nuclei in a detector, per unit detector mass and per unit of time, as a function of the nuclear recoil energy 
E r . The differential event rate per unit detector mass and per unit of time can be written as: 



dN a oPQ a f M f(v x _ 

dE r 2m 2 r m x r J Vtnin ( E ., " X 



F{E r f J -^dv x , (5) 



where <jq is related to the WIMP-nucleon cross-section, a x - p , by <jq = <J x -p ' {Am r /M r ) 2 , with M r = 
" lx , T " p the WIMP-nucleon reduced mass, m r = " lx . mjv the WIMP-nucleus reduced mass, m v the WIMP 
mass, ttin the nucleus mass, and A the atomic weight. F is the nuclear form factor; in the following 
analysis the Woods-Saxon form factor will be used. p ~ 0.385 GeV cm~ 3 and f(v x ) are the density and 
the velocity distribution of WIMPs near the Earth. Let us note that we are assuming identical WIMP- 
proton and WIMP-neutron cross-sections, and that we are ignoring the spin-dependent interactions. In 
our study we will consider a XENON detector with 7 energy bins between 4 and 30 keV. We also consider 
a negligible background. Furthermore, we examine three 'benchmark' experimental setups, assuming 
exposures e — 30, 300 and 3000 kg-year, which could correspond e.g. to a detector with 1 ton of xenon 
and 11 days, 4 months or 3 years of data acquisition respectively. 

Figure [2] shows the exclusion lines (black lines) for exposures e = 30, 300 and 3000 kg-year, on 
the [mo, TO1/2] (upper panels) or the [Mi, fj] (lower panels) parameter space for the two models defined 
previously. The left plots correspond to plain MSSM scenarios whereas the right to the BMSSM, with 
the ei parameter turned on. These curves reflect the XENON sensitivity and represent its ability to 
test and exclude different regions of the parameter space at 95% CL: all points lying below the lines are 
detectable. We note that when some line is absent, this means that the whole parameter space can be 
probed for the corresponding exposure 

As a general rule, the detection prospects are maximised for low values of the pairs m-o and TO1/2 or 
Mi and \x because they give rise to a light LSP. On the other hand, the regions of low TO1/2 or Mi ~ \i 
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Figure 2: Regions in the [mO, mi/2] plane (upper panels) or [Mi, p] plane (lower panels) in which the WMAP constraint 
is fulfilled (red lines). The blue region is excluded by the null searches for charginos at LEP. The orange areas in the upper 
(lower) panels are excluded because the stau (stop) is the LSP. The dash-dotted black curves in the up panel are contour 
lines for the Higgs mass with values in GeV as indicated. Above the dotted purple line of the left of this line, the electroweak 

vacuum is metastable. 



are also preferred because in that case the lightest neutralino is a mixed bino-higgsino state, favouring 
the X? — Xi ~ h and Xi~x1^H couplings, and therefore the scattering cross-section. The detection 
prospects are also maximised for low values of tan /3. The introduction of the NR operators gives rise to an 
important deterioration of the detection prospects. The main effect enters via the important increase in 
the lightest CP-even Higgs mass. Nevertheless, let us emphasize again that this deterioration is relative, 
since we are comparing with a plain MSSM, which is already excluded because of the light Higgs mass. 
Concerning the plots in figure [5J a further remark that can be made is that, even for low exposures, a 
sizable amount of the parameter space can be probed. Larger exposures could be able to explore almost 
the whole parameter space taken into account. 



4 Indirect detection of dark matter with gamma-rays 

The differential flux of gamma-rays generated from DM annihilations and coming from a direction forming 
an angle i/j with respect to the galactic center (GC) is 
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where the discrete sum is over all DM annihilation channels, dN*/dE 1 is the differential gamma-ray yield 
of SM particles into photons, (av) is the total self-annihilation cross-section averaged over its velocity 
distribution, p is the DM density profile, r is the distance from the GC and Bri is the branching ratio 
of annihilation into the z-th final state. The integration is performed along the line of sight from the 
observation point towards the GC. In the following study we take into account three DM density profiles: 
NFW with and without adiabatic compression and the Einasto profile. In our study we will consider the 
FERMI satellite, taking into account a five-year mission run, and an energy range extending from 1 up to 
300 GeV, with 20 logarithmically evenly spaced bins. We focus ourselves to an solid angle Ail = 3 • 10~ 5 
sr around the GC. We will consider two sources for the high-energy gamma-ray background corresponding 
to the galactic emissions coming from resolved and diffuse sources. Both measurements has been done 
by the HESS collaboration. 



In figure [3] we present the detectability regions for the two models described previously, in the 
[mo, TO1/2] (upper panels) and [Mi, /j] (lower panels) parameter space. FERMI will be sensitive to the 
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Figure 3: Regions in the [mO, mi/2] plane (upper panels) or [Mi, p] plane (lower panels) in which the WMAP constraint 
is fulfilled (red lines). The blue region is excluded by the null searches for charginos at LEP. The oranges area in the upper 
(lower) panels are excluded because the stau (stop) is the LSP. The dash-dotted black curves in the up panel are contour 
lines for the Higgs mass with values in GeV as indicated. Above the dotted purple line of the left of this line, the electroweak 

vacuum is metastable. 



regions below the contours and to the areas inside the blobs. The detection prospects are maximised 
for low values of the pairs too and mi/2 or Mi and \i. For the mSUGRA-like model (upper panels), 
the growth of TO1/2 gives rise to the opening of some relevant production channels, after passing some 
thresholds, increasing significantly {uv). The first one corresponds to a light neutralino, with mass 
m x ~ mz/2 ( ttii/2 ~ 130 GeV). In that case the annihilation is done via the s-channel exchange of a real 
Z boson. The second threshold appears for m x ~ mw (mi/2 ~ 220 GeV). The annihilation cross-section 
is enhanced by the opening of the production channel of two real W bosons in the final state. The last 
threshold corresponds to the opening of the channel X1X1 ~~ * tt (toi/2 ~ 400 GeV). The aforementioned 
threshold appears as a particular feature on the upper plots: an isolated detectable region for TO1/2 ~ 400- 
500 GeV and too < 300 GeV. Furthermore, larger values for the annihilation cross-section can be reached 
for higher values of tan/3. For the second benchmark (lower panels), the region where Mi 3> /i is highly 
favored for indirect detection due to the fact that the LSP is higgsino-like, maximising its coupling to 
the Z boson. Let us recall that the Z boson does not couple to a pure gaugino-like neutralino. The 
introduction of the NR operators gives rise to a mild signature when the neutralino is almost a bino-like 
state. However, when /x < Mi there is an important increase of the Xi — Xi — ^ coupling, and therefore 
to a boost in the annihilation into fermion pairs. On the other hand, as the Higgs boson h becomes 
heavier, the processes giving rise to the final state hZ get kinematically closed. Concerning figure [3] let 
us note that the only astrophysical setup in which some useful information can be extracted is the NFW C 
one. This means that in this scenario, in order to have some positive detection in the 7-ray channel, 
there should exist some important enhancement of the signal by some astrophysical mechanism (as the 
adiabatic contraction mechanism invoked in this case). 
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